THE JACOBI MATRICES APPROACH TO NEVANLINNA-PICK PROBLEMS 
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^SJ ' Abstract. A modification of tlie well-known step-by-step process for solving Nevan- 

linna-Pick problems in the class of Rg-functions gives rise to a linear pencil H — \J, 

where H and J are Hermitian tiidiagonal matrices. First, we show that J is a positive 

^ operator Then it is proved that the corresponding Nevanlinna-Pick problem has a unique 

^1^ , solution iff the densely defined symmetric operator J~2HJ~2 is self-adjoint and some 

ciiteria for this operator to be self-adjoint are presented. Finally, by means of the opera- 
tor technique, we obtain that multipoint diagonal Fade approximants to a unique solution 
ip of the Nevanlinna-Pick problem converge to ip locally unifoiTnly in C \ K. The pro- 
posed scheme extends the classical Jacobi matrix approach to moment problems and Pade 
approximation for Rq -functions. 
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ro! . 1. Introduction 

The connection with Jacobi matrices has led to numerous apphcations of spectral tech- 
niques for self-adjoint operators in the theory of moment problems, orthogonal polynomi- 
CN ' als on the real line, and Pade approximation. Let us recall some basic ideas of this interplay. 

First, note that one of the key tools in relating these theories is the class Rq of all functions 
having the representation 



> 



^ : (1.1) ,.(A) . f ""<" 






(<— ^ ■ where cr is a probability measure, that is, J^d<T{t) = 1. If the support supper of <t is 

r^ I contained in [a, /3] we will say that ip S R[a, /3]. 

Consider a probabiUty measure a such that all the moments 



(1.2) s„ := / t"dcr(t), neZ+:=NU{0} 

are finite. In this case, the corresponding function ip has the following asymptotic expan- 
sion 

(L3) <p(A) = --- — __+o|^__j, A^oo, 

for every n e Z_|. (here and throughout in the sequel A^oo means that A tends to oo non- 
tangentially, that is, inside the sector e < arg X < tt — e for some e > 0). In view of the 
Hamburger-Nevanlinna theorem [I], the classical moment problem reads as follows. 
Hamburger moment problem. Is the function ip <E Hq satisfying il.3\t uniquely deter- 
mined by the sequence {sj}°^q of moments? 
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The moment problem is called determinate if Lp is uniquely determined. Otherwise the 
moment problem is said to be indeterminate. In fact, one can give an answer to the question 
in terms of the underlying Jacobi operators generated by Jacobi matrices. To see Jacobi 
matrices in this context, note that one can expand Lp into the following continued fraction 
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where aj are real numbers, hj are positive numbers (see IT], P9l . BOl ). Moreover, numbers 
Cj and hj can be explicitly expressed in terms of the moments sq, . . . , S2j+i ID- Continued 
fractions of the form (11.4b are called J-fractions ll35l . ll49l . To the continued fraction ( 11.4b 
one can associate a Jacobi matrix H and its truncation i?[o,n-i] 
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Let ^?p ^> denote a Hilbert space of complex square summable sequences {xo,xi, 
equipped with the inner product 



ix,y) 



OO 

i=0 



x,yef 



[0,oo)- 



Now, in the standard way, we can define a minimal closed operator H acting in d'^ generated 
by the matrix H IT], IT2l . We will denote the domain of H and the range of H by dom H 
and ran H, respectively. It is easy to see that H is symmetric, i.e. 

{Hx,y) = {x,Hy), x,y Cz doniiJ. 

Moreover, it is well known that H is self-adjoint if and only if the corresponding moment 
problem is determinate and the solution of the problem admits the representation 

^iX)^{{H-X)-'eo,eo) 

where e = (1, 0, ... )^ is a column vector (see HI, 1421 ). In the indeterminate case, a 
description of all ip € Rq satisfying ( 11.31 ) can be found in HI, ifTSl . 1421 (see also l24l 
where the operator approach to truncated moment problems was proposed). In both cases, 
we have 



On (A) 



{{H[Q,n-i] - A) ^eo,eo) 



1 



bl-2 



A - ao 



A — a„_i 



PniX) 

where P„ are orthogonal polynomials with respect to a, and Qn are polynomials of the 
second kind (see HI, l40l . l42l ). It is an elementary fact of the continued fraction theory 
(see, for instance, H], 0, [35]) that 

Qr,(A) „ f 1 



(1.5) 



ifiiX) 



= o 



X- 



PniX) \X^-+\ 

In other words, relation ( 11.51 ) means that the rational function —Qn/Pn is the nth diagonal 
Pade approximant to (^ at oo (for more details on Pade approximants see 15] ). Now, we 
see that in the self-adjoint case, convergence of diagonal Pade approximants appears as 
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the strong resolvent convergence of the finite matrix approximations i?[o,n] to H. So, if 
the moment problem is determinate then the corresponding diagonal Pade approximants 
converge to the solution Lp locally uniformly in C \ R. This statement for the class R[a, /3] 
is known as the Markov theorem BOl . The above-described scheme has been recently 
extended to the case of rational perturbations of Nevanlinna functions ll20l . 11211 . 1221 . 
Also, the scheme was adapted to the case of complex Jacobi matrices ifTOl and generaUzed 
to the case of band matrices fg). 

The main goal of this paper is to generalize the scheme to the case of Nevanlinna- 
Pick problems and to prove convergence of related multipoint diagonal Pade approximants. 
To show our purpose more precisely, let us recall that the classical Hamburger moment 
problem is the limiting case of the following problem (see IT], ll27l . Il36l ). 
Nevanlinna-Pick problem. Let {zfc}^Q be a sequence of distinct numbers from the upper 
half plane C+ and let Lp e Rq. Define numbers Wj := fizj). Is the function cp G Ro 
satisfying the interpolation relation (p{zj) = Wj, j E Z+, uniquely determined by the 

given data {zfej^o- {"'fclfc^o? 

In view of the classical uniqueness theorem for analytic functions, the answer to this 
question is trivial if the sequence {zk}^Q has at least one accumulation point in C+. 
So, in what follows we will suppose that the sequence {zk}1°^Q does not have any 
accumulation point in C+. In other words, all the accumulation points of the sequence 

{zk}T=o lie in «■ 

Similarly to the moment problem case, the Nevanlinna-Pick problem is called determi- 
nate, if ip is uniquely determined. Otherwise the Nevanlinna-Pick problem is said to be 
indeterminate. We should also note that diagonal Pade approximants at oo are the limiting 
case of the following multipoint diagonal Pade approximants. 

Definition 1.1 (Q). The nth multipoint diagonal Pade approximant for the function Lp at 
the points {zq, Zq, . . . , z-,-, z-,-, . . . } is defined as a ratio —Qn/Pn of two polynomials Qn, 
Pn of degree at most n — 1 and n, respectively, such that the function Pnf + Qn vanishes 
at the points zq, 'zq, ■ ■ ■ , Zn-i, Zn-i- 

It appears that the problem of finding multipoint diagonal Pade approximants for the 
Ro-function ip at the points {zq, zq , ■ ■ ■ , zjjJj, . . .} is closely related to a continued frac- 
tion expansion of the following type 

(l^g^ 1 I bliX-zo){X-zo)\ 62(A-zi)(A-zi)| 






a?h-a^^ 



4^)A-a« 



where a, are real numbers and a^ , bj are positive numbers. This continued fraction 
gives rise to a tridiagonal linear pencil H — XJ, where H and J are semi-infinite tridiag- 



onal matrices 11231 (see also 11501 where tridiagonal linear pencils associated with general 
continued fractions of type ( 11.61 ) were introduced). In this paper, we firstly obtain that J 
generates a positive operator Then we introduce a densely defined symmetric operator 
J^2 HJ^2 and present criteria for this operator to be self-adjoint. Next, we prove that 
the NevannUna-Pick problem in question has a unique solution if and only if J^^HJ~2 
is self-adjoint. Finally, we show that if J^^HJ^^ is self-adjoint then the locally uniform 
convergence of the multipoint diagonal Pade approximants 

Qri+l(A) _ /, -i -i , 1 -i -i N 
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to the unique solution 

y>(A) = ({J---HJ-"^ - x)-'J-ho, J-ho) 

of the Nevanlinna-Pick problem arises as the resolvent convergence. 

The paper is organized as follows. In Section 2 we present the step-by-step process 
for solving the Nevanlinna-Pick problems and associated sequences of polynomials. In 
Section 3, a tridiagonal linear pencil is introduced and basic properties of the operator 
J are given. The one-to-one correspondence between tridiagonal linear pencils and the 
Nevanlinna-Pick problems in question is shown in Section 4. The next session is concerned 
with the Weyl circles. Section 6 reveals the underlying symmetric operators. In Section 7, 
we characterize the determinacy of the underlying Nevanlinna-Pick problems in terms of 
the self-adjointness of J^^HJ^^. After that, in Section 8, for the determinate case, we 
prove the locally uniform convergence of multipoint diagonal Pade approximants for Ro- 
functions. 

2. The modified multipoint Schur algorithm 

As is known, the Schur transformation is a powerful tool in solving moment and in- 
terpolation problems (see |[T|, lO). The starting point of our analysis is the following 
modification of the Schur transformation. 

Proposition 2.1 (cf. II23I ). Let ip G Ro and let z G C+ be a fixed number. Then there exist 
unique numbers a^^', a^'^' € M and b > such that the function ipi defined by the equality 

^^■^^ ^^^^ " " a(2)A - a(i) + 62(A - z){X - z)(^i(A) 

belongs to Rg U {0}, that is, ipi has the representation (II. Il l with a probability measure in 
case ipi ^ 0. Moreover, we have that 

(2.2) b^ = a(2) - 1. 

Proof. To see that the numbers a'^', a'^^ are uniquely determined, let us substitute A for z 
and z in ( 12. 11 1. We thus get 

(2.3) a(2)^_„(i) = _ 1 , a(2)--aW- — 



ip{z) Lp{z) 

Eliminating from the above relations a*^^' and a'-^^ one can obtain the following formulas 



(2.4) a 



(1) 



tda{t) 



da{t) 



a(2) 



d(7{t) 



da{t) 



\t-z\\ 
Further, it follows from the Schwartz lemma that 

{\-z)(\-z) Jmt-X 

(the proof of this fact is in line with that of |[23l Lemma 3.1]). Choosing 6 > in the 
following way 

b^ = f dn{t) 



and defining (pi := (fi/b'^ we get that the function ipi possesses the integral representa- 
tion dl.lt with a probability measure. Finally, by taking X ^ iy and y — ^ oo in ( |2.5l l we 
get (O. D 
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Remark 2.2. It should be noted that for Lp e R[a, /3] this modification of the Schur algo- 
rithm was presented in 1231 Lemma 3.1]. However, its proof is valid for Lp £ Rq. A similar 
transformation for Caratheodory functions was proposed in |fT9l . 

Let (/9 be a non-rational function of the class Rq, i.e. (p admits the representation ( II. Il l 
with a probability measure which has an infinite support. Let also an infinite sequence 
{^k}kLo "^ '^+ of distinct numbers be given. Since ip is not rational the given data give 
rise to infinitely many steps of the step-by-step process. So, we have infinitely many linear 
fractional transformations of the form ( 12.11 ) which lead to the following continued fraction 

(2.6) 1 I bl{X-zo)i\-zo)\ 52(A-zi)(A-zi)| 






af'A-aW 



a^^h-ai'^ 



(for more details, see li23l ). It should be noted that general continued fractions asso- 
ciated with finding multipoint Pade approximants were introduced in ||32|| and studied 
in I33l, m. 

It is immediate from the construction that the {n + l)th convergent of ( 12.6b 

l(A) 1 I 6Li(A"^«-i)(A-z„-i)| 



^n+U 



Pn+l{X) 
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satisfies the following interpolation relation 

(2.7) ^(zj) = -- — -, j^O,...,n. 

Since (p € Rq and the coefficients a^ ' , a^- ' , bj are real, one also has 

/_ N Qn+l{Zj) . „ 

^{zj) = —fS 7=^^ j=0,...,n. 

^7l+l[Zj) 

So, we have just concluded the following. 

Proposition 2.3. The rational function ~Q„+i/Pn+i is the (n+l)th multipoint diagonal 
Fade approximant to ip at the points {zo,'zo, ■ ■ ■ , Zj^'Zj, . . .}. 

It is well known that denominators and numerators of convergents of a continued frac- 
tion satisfy a three-term recurrence relation (see, for instance, ||35]| ). In particular, for the 
continued fraction ( 12.61 ) the recurrence relation takes the following form 

(2.8) uj+i ~ [afx - af)u, + b]_^{\ - 2j-i)(A - z,-i)m,-i - 0, j G N. 

Further, the polynomials Pj of the first kind are solutions Uj = PjW of the system ( 12. 8t 
with the initial conditions 

(2.9) Mo = l, ui = 4'U-a|,'^ 

Similarly, the polynomials of the second kind Qj{X) are solutions Uj = Qj{X) of the 
system ( 12. 8t subject to the following initial conditions 

(2.10) Mo = 0, ui = -1. 



Remark 2.4. Note that the polynomials Pj are orthogonal with respect to the varying 

dcr(t) 

an operator treatment of the relation of the polynomials Pj to orthogonal rational func- 
tions was presented in ||23l (see jfTTl Section 9.5], where this relation is also discussed). 



measures — ..^ ^ ' (see 1291, E), gTl Section 6.1]). Moreover, for (p e R[a,^] 
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It should be also remarked that some orthogonality relations for polynomials and rational 
functions related to general continued fractions of type ( 12.61 ) were obtained in 1341 . lISOl 
(see also [jSTTl, where biorthogonality properties of rational functions related to multipoint 
Pade approximation were studied and concrete examples connected with generahzed hy- 
pergeometric functions were constructed). 



3. Tridiagonal linear pencils associated with Ro-FUNCTIONS 

In order to see linear pencils in our context, let us note that the recurrence relation ( 12.81 1 
can be renormalized to the following one 



(3.1) 



(bj_i - Ac)j_i)mj„i + (uj - Xcj)uj + {bj - Xdj)uj+i = 0, j e N, 



where the numbers Oj, bj, Cj, Oj are defined as follows 



ij =«f\ 



""A 



J' 



(2) 



dj ~ bj , j e z_| 



and the transformation w — > u has the following form 

^ ^ Ui 



(3.2) 



Uq = Uq, Uj 



bo... bj^i{zQ - X) ... (zj_i - A) ' 



jeN. 



Thus, we have two associated sequences Pj and Qj of rational functions obtained from the 
polynomial sequences Pj and Qj, respectively, by means of the transformation (13.2b . In 
contrast to the polynomial case, the rational functions Pj are not orthogonal with respect 
to the original measure a since 



(3.3) / Poit)Pi{t)dait) = [ Pi{t)da{t) 



J2) 



and, due to ( IZ4l i. 1 



,(2) 



7^ for any zq £ C+. Despite this, some orthogonality prop- 



erties remain valid (see ifTTl Theorem 2.10]). It should be also noted that some orthogonal 
proper rational functions satisfy a relation similar to (13.11 1 lH p. 541] (see also ifTTl for the 
recurrence relations for orthogonal rational functions). 

The relation ( 13.1b naturally leads to a linear pencil H — XJ, where 



/ Op bo 
bo ai bi 



H = 



02 



V 



\ 



/ 



J 



/CO 

5o 



V 



Cl 



5i 

C2 



\ 



/ 



are Jacobi matrices. For an infinite matrix A, we denote by A[j ^.j the square sub-matrix 



obtained by taking rows and columns / 
and k we have that 



j,j + l,...,k < oo. For example, for finite j 
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flfcy 



-^UM - 
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By J we also denote the minimal closed operator on £?p ^^ generated by the matrix J |[T]. 
Obviously, J is a symmetric operator Besides, due to ( 12. 2b . we have the relation Cj = 
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1 + d^, which gives us the following factorization of J 
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J ^ L*L 
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The factorization of J allows us to say a bit more about J. 

Proposition 3.1. TVie operator J is self-adjoint and positive, that is, 

{Jx,x)>0, xedomJ\{0}. 

In particular, ker J = {0}. 

Proof. Let us consider the Hermitian form ( Jf , ^) on finitely supported sequences ^, that 
is, ^ = (^0,6, ■ • ■ iCmOjO, . . .)^. By virtue of (13.411 . we have that 

Further, let us prove that ker J* = {0}. Suppose the converse, that is, there exists 7] e i'^ 
such that J*ri = and ?/ 7^ 0. Taking into account the structure of J we get the equality 

= (J*?7, 77) = |?7oP + lOo* + ?7ip + . • • + |5„_i77„_i + 77„|2 + . . . ^ 

which implies 77 = 0. So, ker J = ker J* = {0}. This contradiction also shows that 



(3.5) 



CXD 

E 

fe=0 



bfc 



where pj are polynomials of the first kind associated with J. Since the relation ( 13.5b 
doesn't hold true for Jacobi operators with deficiency indices (1,1) (see IIT2II . ll42ll ). we 
obtain that J is self-adjoint. The statement of the proposition also immediately follows 
from mi Theorem VII. 1 .4] . D 

Remark 3.2. It has been recently proved ifTTl that if (p G R[q;, /3] and Zk —)' 00 then 

{Jx,x) > S{x,x), X e £^, 

for some 6 > 0. Furthermore, in this case the operator J is a compact perturbation of / and, 
in fact, the linear pencil H — XJ is a compact perturbation of the classical pencil Hq — XI 
(which corresponds to the limiting case Zk ^ 00 for A; = 0, 1, 2, . . . ). It should be noted 
that in the case of orthogonal Laurent polynomials a similar tridiagonal pencil was consid- 
ered in ITSl . Roughly speaking, the case of orthogonal Laurent polynomials corresponds 
to the multiple interpolation at and 00, which is known as the strong moment problem on 
the real line |[35l . An operator approach to the strong moment problem was given in ||3TI . 
It is also worth to note that, in the matrix case, Jacobi type symmetric operators related to 
the matrix strong moment problems were presented and studied in ll44l . ||45]| . 

Since ker J = {0} and J is self-adjoint, we can consider the self-adjoint operator J~ 2 , 
which is not necessarily bounded. However, the following statement holds true. 



Proposition 3.3. We have that 
(3.6) e, 

where the vectors eo = (1, 0, 0, . 

in e. 



G dom J 



(0,1,0,...)^ 



form the standard basis 
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Proof. It is the basic spectral theory that for the positive operator J there exists a resolution 
of the identity Et such that 

/•oo 

J/= / tdEtf, /edomJ, 
Jo 

and / S dom J if and only if j^ t'^d{Etf, /) < oo 111 Section 66]. Moreover, we also 
have that 

1"°° 1 
J-^f^j -FdEtf, fedomJ-^, 

and / G dom J~ 2 if and only if L jd{Etf, /) < 00. Now, (13.61 ) is equivalent to 

.00 ^ 

/ -d{Etej,ej) < 00, j e Z+. 
Jo t 

First we will prove that 

r°° 1 

(3.7) / --d{Eteo,eo) < 00. 

Jo ^ 

For simplicity, let us denote ly = {E.eo,eo) and introduce the similar measures i/„ = 

(£^.^"■^60, eo) for the truncations J[o,n], where i? is such that 

/>oo 

-^[o.n] = / tdE^''\ nel.+ . 
Jo 

Next, it is a standard fact of theory of moment problems [I], that 

/•oo 

'ip{t)dvn{t) ^ / ip{t)dv{t), n-^oo, 
Jo 

for any simple function ip (that is, tp is measurable and assumes only a finite number of 
values). Now, recall that in ll23l Lemma 6. 1] it was proved that 

(3.8) / ^di^nit) = (^[o,l]eo,eo) < 1, ne Z+. 

Thus, Fatou's lemma for varying measures BTl Proposition 17, p. 231] and ( 13.81 ) yield 

(3.9) / -diy{t) < liminf / -di/„(i) < 1. 

Jo t n^oo _/p t 

The rest is a consequence of ( |3.7| i. Indeed, it is well known that for any A from the resolvent 
set p{J) of the operator J we have the following formula for the diagonal Green function 

(3.10) {{J-X)-'e„e,) ^ p,{X) {pj{\) {{J - X)-'eo,eo) + q,{X)) , j € Z+, 

where pj and qj are polynomials of the first and second kinds, respectively, associated with 
the Jacobi matrix J (see for example ifTOl Theorem 2.10], 1281 Proposition 2.2]). Putting 
A = —X, X > 0, into formula ( 13.101 ). it can be rewritten as follows 

-d{Etej,ej) = Pj{-x) I Pj{-x) / — — diy{t) + qj{-x) ) , j G N, 



Q t + X y Jq t -\- X 

where pj{—x) ~ ° |^' '°'iT^' — — > for a; > 0. Now, it remains to apply the Fatou lemma 
to Jq°° j^d{Etej , Cj) as X ^ and to use (13.7b . D 
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Remark 3.4. The main ingredient in the proof was to obtain (13.7b . Another way to prove it 
is through the Darboux transformations. Namely, let us consider a Jacobi matrix Ji = LL* 
and let v* be a corresponding probability measure associated with Ji. Then it follows 
from llT6l Theorem 3.4] that 

du{t) = ctdiy*{t), c> 0. 

The latter relation immediately implies ( 13.71 ). 

To end this section, note that we can now say more about the sequence I >/ro^ieoj eo 
Namely, the following relation holds true 

(3.11) (^Jjp_^„jeo,eoj -> 1, as n^oo. 

Indeed, by applying flE', formula (2.15)] we see that I J,^^^,eo, eoj, n £ Z-|_, are conver- 
gents to the continued fraction 

1| OqI 5i| 

r^ fcT fcT 

Forasmuch as Cj = 1 + 5|, applying the remark to Sleszyiiski-Pringsheim's theorem given 
on 135] p. 93] implies dim . 

4. Relations between Nevanlinna-Pick problems and linear pencils 

In this section we show that there exists a one-to-one correspondence between the linear 
pencils under consideration and the Nevanlinna-Pick problems in question. We also re- 
examine some facts for the polynomials Pj and Qj which are well known for orthogonal 
polynomials. 

We begin with the following connection between the polynomials of the first and second 
kinds Pj, Qj and the truncated linear pencils XJ[qj] — ^[Oj]- which in the classical case 
can be found in 1121 Section 7.1.2] and ||4] Section 6.1]. 

Proposition 4.1. The polynomials Pj and Qj, j G N, can be found by die formulas 

(4.1) Pj{X) = dct(AJ[oj_i] - H[Q,j_^), Qj(A) - dct(AJ[i,,_i] - H[i,,_i]). 

The zeros of the polynomials Pj and Qj are real. Moreover, the polynomials Pj and Qj 
do not have common zeros. 

Proof. Formula ( 14.11 ) immediately follows from the definition of Pj and Qj by using the 
Laplace expansions of the determinants by the last row. Since J[oj-i] is strictly positive, 
one can rewrite the first relation in (14. It as follows 



Pj{X) = det J['o(^'_^] dct(A - J[oy_'i]^[Oj-i]^[oy-i]) det J[oj_i]. 

Clearly, J,g _-|^, iJ[o,j_i] Jrg _ j^i is a self-adjomt matrix. Thus, the latter relation yields the 
fact that the zeros of Pj are real. Similarly, one can show that the zeros of Qj are real. The 
last statement follows by induction via applying the Laplace expansion of the determinant 

det(A J[o j-i] — i?[o j-i] ) by the first row. D 

By induction, one easily gets from (13.1b the Liouville-Ostrogradsky formula 

n-l 

(4.2) Q„+i(A)P„(A)-Q„(A)P„+i(A) = [] 6^(A - Zfc)(A - z^), 

fe=0 
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for every n E Z+ (see ifTTI '). Going further in this direction, we should note that, some- 
times, it is very useful to have (13. Il l in the following matrix form 

(4.3) {H - AJ)7r[o,,](A) = -(b, - Ac),)P,+i(A)e, + (fa, - AD,)P,(A)e,+i, 

(4.4) (H - XJ)C[o,j] (A) - -(fa, - At),)4+i(A)e, + (fa, - A5,)Q,(A)e,+i + eo, 
where the vectors ttjo,] (A) and ^[q ,] (A) are defined as follows 

7r[o,,](A) = (Po(A),Pi(A), . . . ,P,(A),0,0,. . .)^ , 



e[o.j](A) = (Qo(A),Qi(A),...,Q,(A),0,0,. 



For example, by virtue of ( 14.31 ) we get the following generaUzation of the Christoffel- 
Darboux formula. 

Proposition 4.2. We have that for j G Z+ 



(A-C)5](Pfc(A) + c)fc_iPfc_i(A))(Pfc(C)+5fc_iP,._i(C))- 

(4.5) ''=° 

_ P,-+i(A)P,-(C)-P,+^(C)P,(A) 

where 0-1= Ofor convenience and A, C G C-|- \ {zfc}|,^, 



fc=o- 



Proof. It clearly follows from (14.3b that 



(4.6) {{H - AJ)7r[o,,](A),^[o,,](C)) = -(fa, - AO,)^,+i(A)P,(C), 



(4.7) ((i7-a)^[o,,](A),^[o,,](C)) = -(fa, -C5,)^.-+i(C)^,(A). 

Subtracting jA-M from ( 14. 7t and using ( 13.21 ) we get the following relation 

(4.) ,A-o(,/.,M(A),.„(0) ^ ^n'-vg'";fT - 

nA;=o^fc(A-^fc)(c-2fe) 

Now, observe that due to (13.4b we have 

{jn[o^j]{X),TT[oj]iO) ^ (P7r[oj](A),i7r[oj](C)) 
and, so, from (14. 8b we obtain ( 14.51 ). D 

Remark 4.3. To see how it is related to the classical Christoffel-Darboux relation IT] 
let us note that, according to ( 12.4b and (12.2b . we have that dk -^ and bf./\zk\'^ — > 
b1 y^ as Zk —>■ oo, k = 0,. . . ,j provided that the numbers J^t''da{t) are finite for 
k = 0, . . . ,j. Consequently, the classical Christoffel-Darboux formula is the limiting case 
of ( 14. 5b . Moreover, it is shown in 1231 Theorem 2.2] (see also ifTTl Section 4]) that the se- 
quence {Pfc+5fe-iPfc-i}fc^o '^ ^ sequence of rational functions orthogonal with respect to 
the original measure a (see ifTTl for further information on orthogonal rational functions). 
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In what follows we will also need the following relation 

3 
k=0 

(4.9) = {J{^^%j] (A) + ^[oj] (A)), (w^[o,,] (A) + C[o,,] (A))) - ' 



A- A 



1 \ojP,{\) + Q,{\)\\^^^ujPj+^{X) + Qj+i{\) 



Ii^ ^ nCo ^fe I A - zk 1 2 """ i^P, ( A) + Q, ( A) 

where lo G C+ and A e C+ \ {zk}l.^Q- Formula ( |4.9t can be easily obtained by straight- 
forward manipulations with ( 14. 3t and ( 14.41 ) (for the classical case see IT] Section 1.2.1]). 
Next, by following ll28l . let us introduce 7Ti-functions of the truncated linear pencils. 



Definition 4.4. Let j and n be nonnegative integers such that j < n. The function 

(4.10) ™b-,Ti](A) = ((ff[j,,i] - XJ[j^„])~^ej,ej) 

will be called the m-function of the linear pencil H[j,n] ~ ^J[j,n\- 

To see the correctness of the above given definition it is sufficient to recall that J[j.„] is 
positive definite in view of Proposition l3. 1 l and to rewrite ( 14.10b in the following form 

(4.11) m[,.„](A) = ((J^!ji/[,,„] J[^i] - A)-VfTjje„ J'^^e^ . 



Literally as in the classical case (see for instance 112811 ). one obtains that m-functions satisfy 
the Riccati equation. 



Proposition 4.5 ( II23II ). The m-functions rnijn] ^'^d r/irj+i^n] ^^^ related by the equality 

1 



,(2)\ Jl) 



(4.12) '^[j,n] = 

Oj^^'A - ay + bj{X - Zj){X ~ ZJ)"^[J+l,„] (A) 

The latter statement allows us to see the relation of m-functions to multipoint diagonal 
Pade approximants. 

Proposition 4.6. Let dn = dot J[o,,i]/dct J[i „] and r]n = dct J[o,n]/dct J[o,ii-il- Then 
the function Onm^Q^n] '■* '^n Ro-function and 

,A T7\ /\\ Qn+l(A) 

(4-13) "^[o,«](A) = -^ TT^, 

that is, 77i[o.n] is the (n + l)th multipoint diagonal Fade approximant for if. Moreover, we 
have that —rjnPn / Pn^i G Rq. 

Proof. Formula (14.13b is implied by the relation (14.12b . Now, from Proposition 12.31 we 
see that ?7i[o,n] is the [n + l)th multipoint diagonal Pade approximant for ip. To see that 
^n?^i[o,n] G Ro, it is enough to recall that <1> e Rq if and only if 

Im$(A) 
ImA 
and sup [^^(iy)! = 1 [H Section III. 1.1]. The first condition is easily verified by means 

of ( 14.111 ) and the second one follows from ( 14.131 ). In the same way, by noticing that 

Pn{X) dct(AJ[n,„_i] -i7[o,„-i]) , X 7 N-i \ 

~~5 7T\ = 1 ,t\ 7 u ^ — = U-W[o.n] - AJ[o,„] e„,e„j 

Pn+l(.X) dct(AJ[o,„] -H[o,„]) 

one can check that —rjnPn/Pn+i G Ro since 7^„ > 0. D 
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Due to -9nQn+i/Pn+i G Ro and -rjnPn/Pn+i e Ro, we get the following. 
Corollary 4.7. We have that 

i) The zeros of Qn+i and Pn+i interlace, 
ii) The zeros of Pn and P„+i interlace. 

Summing up Propositions l2. Il and l4.61 we conclude the following. 

Theorem 4.8. There is a one-to-one correspondence between the linear pencils in question 
and the data {zk}'^Q, {wk}'^Q of the Nevanlinna-Pick problems. 

Proof. It follows from formulas ( 12.21 ) and ( 12.41 ) that the data {zfcj^g, {wfej^o uniquely 
determine the linear pencil, that is, the following numbers 

(4.14) aj ^ af\ bj ^ Zjbj, Cj = a^^ , d^=bj, j G Z+, 

where a^- € IR, a^- > 0, bj > 0, Zj € C+, and Cj = 1 + 5^. Let us suppose that we 
are given a set of numbers that can be represented as above. Then we see from ( |4.14| i that 
Zj = bj/()j. Finally, by virtue of Proposition |4.6| we get that the numbers Wj are uniquely 
determined by the formula 

^^" Pniz,) 

for large enough n. It remains to note that in view of the precompactness of the family 
—Qn/Pn (see Proposition 18. lb and ( 13.11b there exists a function Lp E Rq which satisfies 
the underlying interpolation relation (p{zj) = Wj, j E Z+. D 

5. The Weyl circles 

The classical Weyl circles approach to Nevanlinna-Pick problems can be found in 1271 
Section IV.6]. In this section, following fl\ Section 1.2.3], we adapt the notion of the Weyl 
circles to the linear pencil case. 

Let us begin by considering the function 

where AeC\M, reMU {oo}, and j € N. Obviously, from the definition we have that 

U!j{X,Oo) = U!j-l{\,0). 

Moreover, in view of ( 12.7b we have that Ldj{zk,T) — w^ and Ldj{'zk,T) — Wk for j = 
A: + 2, fc + 3, . . . . So, formula ( 15.11 ) gives a parametrization of [j-l/j] rational solutions to 
the truncated Nevanlinna-Pick problems. Another such a parametrization is given in ifTTl 
Theorem 6. 1 .3] in terms of orthogonal rational functions of the first and second kinds. 

Due to Proposition |4.6l the number — 'p]\\ is not real for any A £ C\M and, therefore, 
we see that the set 

Kj{\) = {ujj{\,t) : T e MU {oo}}. 

is a circle. In addition, we have that Kj{X) = Kj{\). So, we can consider only the case 
when A e C+. The following statement contains a characterization of the circle Kj{X). 

Theorem 5.1. Let A e C+ \ {zkYi^i^ be a fixed number Then the center of Kj (A) is 



.5 2) _ Q,(A)P,--i(A)-Q,-_i(A)P,(A) 

P,(A)P,_i(A)-P,_i(A)P,(A)^ 
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and the radius of Kj (A) is 
(5.3) ^ 



Besides, the equation of Kj{X) can be represented as follows (setting c)_i = Oj 

j-i __ 

(5.4) Y. \^iPk{\) + dk-iPk-iiX)) + Ofe(A) + c)fe-iQfc-i(A)|2 - ^^ = 0. 
fc=o A - A 

Proof. By the same reasoning as in the proof of HI Theorem 1.2.3] we conclude that 

Q,(A)P,_i(A)-gj_i(A)P,(A) 



'^j(A,t) = -■ 



(A)P,_i(A)-g,.i(A)P,(A) 



P,(A)P,_i(A)-P,_i(A)P,(A) 



P,(A)P,_i(A)-P,_i(A)P,(A) 

where 9 ~ 9{t) is real. The latter relation immediately gives us ( 15.21 ) and the formula for 
the radius of Kj (A) 

Q,(A)P,_i(A)-Q,_i(A)P,(A) 



P,(A)P,_i(A)-P,_i(A)P,(A) 

which by means of ( 14.2b and ( 14.51 ) can be reduced to ( 15.31 ). 

The rest of the proof is identical to the proof of (1, Theorem 1.2.3]. D 

Denote by Kj(A) the closure of the interior of Kj{X). Then the following statement 
holds true. 

Corollary 5.2. Let A e C-|- \ {^A;}^Zq be a fixed number Then the set Kj(A) is a set of 
numbers w G C satisfying the inequality 

j-i __ 

(5.5) Y. \^iPkW + ^k-iPk-i[X)) + Qfe(A) + c)fe_igfc_i(A)|2 < ^^^. 

fc=o ^ " ^ 

Furthermore, we can get a relation between the discs K^+i (A) and Kj (A). 
Corollary 5.3. We have that 

K,+i(A)CK,(A), jeN. 

Besides, the circles Kj+i (A) and Kj (A) have at least one common point. 

Proof. The proof of the both corollaries is in line with the proof of the analogous state- 
ments given in IT] Section 2.3]. D 

Now, we see that there are two options for the sequence Kj(A). Namely, we can have a 
limit point or a limit circle. 

Theorem 5.4. Let A e C+ \ {zfej^o ^^ a fixed number Then the sequence Kj(A) 
converges to a point iff 



^|Pfe(A) + i)fe_iPfe_i(A)|2=oo. 



fc=0 



Proof. The proof is immediate from Corollary 15. 3 1 and ( 15.31 ). D 

Next, we obtain the existence of the Weyl solution. 
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Theorem 5.5. For every A G C+ \ {^fej^g there exists a number to = a;(A) G C+ such 
that 

oo 

(5.6) Y. \^iPkW + ^k-iPk-iW) + QkW + ^k-iQk-iW\^ < 



^-^ 



Proof. The statement is a straightforward consequence of Corollary 15.31 and the inequal- 
ity (E). □ 

Finally, it should be noticed that the mentioned parametrization from ifTTl leads to a 
slightly different but very similar theory of nested disks ifTTl Section 10]. That theory is 
equivalent to the presented one in the sense that the underlying Nevanlinna-Pick Problems 
are the same. 



6. The underlying symmetric operators 

In this section we reduce the linear pencil in question to an operator generated by the 
formal matrix expression J^2iyj^2. Namely, we show that this operator is a densely 
defined symmetric operator 

Since ej G dom J C dom Ja the vectors fj := J'ej, j e Z+, belong to i'^. The 
relation ker J 2 = {0} implies that the linear span 

J" = span{/j}°^o ^ 1 ^ ^''■^'' : Cfc e C, n e Z+ l 
I fe=o J 

is dense in d'^. In view of ( I3.6l l. we can also introduce the vectors gj := J^^ej, j G Z+, 
which lie in £^. Moreover, the linear span Q = spa.n{gj}°^Q is dense in i'^. Besides, we 
have that that the systems {fj}°^Q and {gj}°^Q are bi-orthogonal, i.e. 

As a consequence, we get that there is a one-to-one correspondence between h ^ P and 

the formal series 

00 00 

Y{h,gk)fk, Y{h,fk)gk- 

fc=0 k=0 

In this case, we will write h ^ J2kLni^'9k)fk or h ^ X^fcLoC*' fk)9k- Next, we see that 
(setting b_i = for convenience) 

J-'sHJ-^f^ = bj-igj-i + a-jgj + bj.gj+i, j G Z+. 

So, we have that J^^HJ^^ : T 1-^ Q. Thus the domain of the matrix expression 
J^^HJ^2 is dense in £^. 

Proposition 6.1. The formal matrix expression J^^HJ^^ generates a densely defined 
symmetric operator with the deficiency indices either (1,1) or (0,0). 

Proof. It is easy to see that 

{J-^HJ-ifjJk) = {f,,J-^HJ-ifk), J, fc G Z+, 
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that is, J^2HJ^2 is symmetric in P. Thus, the operator is closable and, in what fol- 
lows, by J^^ HJ^2 we denote the minimal closed operator defined by the matrix expres- 
sion J~2HJ~2. Let {J~^HJ~2y be adjoint to J~^HJ~2 in i^. By the definition, a 
vector h G dom( J~ 2 HJ~ 2 )* if and only if there exists a vector h* G £'^ such that 

{j~-2HJ~-^fk,h) = (/fc, h*), fekez+. 

Further, it can be rewritten as follows 

[bk-igk-i +afc5fc + bkgk+i,h) = {fk,h*), k e Z+, 
which actually implies that 

yk = bk~iXk-i + akXk + bkXk+i, k e Z+, 

where h ^ Y^^=o^kfk and h* ^ J2T=oyk9k- Thus, h e dom(J^2iJJ^2 )* if and only 
if there exists h* £ P such that 



h* ^ '^[bk-iXk-i + akXk + bkXk+i)gk- 



k=0 

The next step is to determine the deficiency indices. In order to do that we should find 
nontrivial solutions of the equation 

(6.1) ((J-5HJ-5)*-A)/i = 0, ImAT^O. 

Let h ^ X^fcLo ^kfk be a solution to ( 16.1b . Then we obviously have that 

ifk,iiJ'^HJ-iy-X)h)=0, fceZ+, 
which reduces to the following 

bk-iXk-i + akXk + bkXk+i = X{fk,h), fc e Z+. 

Observing that (A-, h) = dk-iXk^i + tkXk + ^kXk+i, we arrive at 

{bk-i - Xdk-i)xk~i + {ak - \tk)xk + {bk - Mk+i)xk+i =0, fee Z+. 

In view of ( 13.11 ), ( 13. 2t , and ( 12.91 ), we conclude that xk = cPfc(A). So, the linear space M\ 
of the solutions to (16.1b has dimension 1 if there exists an element h € P such that 



(6.2) /i-$]Pfe(A)/fc. 



fc=0 



Otherwise, the linear space M\ has dimension 0. 

Let us find the condition for h from ( 16.21 ) to belong to l"^. First, we should check 

the weak convergence of the sequence h^ ~ X^I—o PkWfk- Obviously, we have that 

{hn,gk) -> {h,gk) = PkW as ri -^ 00. Furthermore, Q = span{,gj}°^Q = £'^. Con- 
sequently, according to the criterion of the weak convergence we get that the convergence 
of ( 16. 2b is implied by the uniform boundedness of the following sequence 



Y.P'^WfkW = (J^[o,„](A),^[o,„](A)) = 

n 

(L^[o,„] (A), Ln^o.n] (A)) = ^ |P, (A) + dk-iPk-i (A)|^ 



(6.3) '"~° 



k=0 
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From (16.31) we see that the condition 

oo 

(6.4) ^|Pfe(A) + c)fe_iPfe_i(A)|2<oo 

k=0 

guarantees the existence of h satisfying (16.21 ). It turns out that this condition is also nec- 
essary. Indeed, let us suppose the converse that J2T=o l^kW + i'/£-i-Ffe-i(A)P = oo 
and there exists h € P having the representation ( 16.21 ). Then it follows from ( 16.11 ) that 
h E ran J^2 and, therefore, h ~ J^ho for some Hq € £'^. The latter means that 

oo 

\\h\\ = WjhioW = \\Lho\\ = Yl \PkW + 0fc~ift-i(A)|2 = oo, 

k=0 

which yields the contradiction. So, dim A/a = 1 if and only if ( I6.4l i holds true. 

It is well known that for symmetric operators the deficiency index d\ = dim A/a is the 
same for each A £ C^ as well as for each A € C_. Further, it follows from ( 15.31 ) that 

n — 1 n — 1 

Y, \PkW+Vk-iPk-iW\^ - Y i^fcW + ^k-iPk-iiW 

k=0 k=0 

since the radii of Kn (A) and Kn (A) are equal. The latter relation implies that d\ = dj. D 

Now we are in a position to formulate criteria for J^^HJ^^ to be self-adjoint (for the 
classical case see |[T], llT2l . B2l ). 

Theorem 6.2. The following statements are equivalent: 
i) The operator J~^ HJ~^ is self-adjoint; 

ii) The sequence Kj(A) converges to a point for some A G C^ \ {zk}^^^; 
Hi) We have that 

oo 

(6.5) ^|Pfe(A) + c)fc_iPfc_i(A)|2 = oo 

k=0 

for some A G C+ \ {^fcj^o- 

Proof. The equivalence of ii) and iii) is established in Theorem 15 .41 The equivalence of 
i) and iii) is actually proved in the proof of Proposition 16.11 by showing that the defect 
vector ( 16.21 ) belongs to i"^ if and only if ( |6.5l l holds true. D 

Remark 6.3. It is well known that for symmetric operators the dimension of the defect 
space A/a remains the same for all A G C+ . Thus, if ( 16.5b holds for some Ao G C+ \ 
{zfej^g then it holds for all A G C+ \ {zfe}^Q. The same is true for the limit point case. 

We should emphasize that in our approach the operator J^^HJ^^ plays exactly the 
same role as the Jacobi matrix for a moment problem. We should also stress here that if 
the original measure has finite moments of all nonnegative orders and we have a collection 
of interpolation sequences {z^ s'kLo such that for every k G Z+ 

Zj, — )■ oo, as n — > oo, 

then the corresponding matrices J*^") converge to the identity /, as n — > oo, elementwise 
(see ( 12.21 ) and ( I2.4l i). So, roughly speaking, in this case, the operator ( j'") ) ^ 2 ij(") ( j(") ) ^ 2 
approaches the classical Jacobi matrix (see also ifTTIl '). 

To complete this section, it should be remarked that, in recent years, a lot of attention 
has been paid to the study of orthogonal polynomials on the unit circle via the spectral 
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theory of CMV-matrices (see H43II and references therein). Roughly speaking, orthogonal 
polynomials on the unit circle correspond to the multiple interpolation problem at and cx) 
for the Schur class (actually, there is only one interpolation point since oo is symmetric to 
with respect to the unit circle). The multiple interpolation at two points is, in some sense, 
the limiting case of the case under consideration. Also note that an operator approach to 
orthogonal rational functions on the unit circle via CMV matrices can be found in BHl . It 
is also worth mentioning that Jacobi type normal matrices associated to complex moment 
problems were introduced and studied in ifTJl . lfT4l . 

7. The uniqueness of Nevanlinna-Pick problems 



In this section, by mimicking the proofs of B2l Theorem 2.10] and 142] , Theorem 2.11], 
we characterize the determinacy of the Nevanlinna-Pick problems in question in terms of 
the self-adjointness of J^^HJ^^, 

Let ip E Rq and let a sequence of distinct numbers {^tj^Q C C+ be given. According 
to ( 12.3b and (12.2b . the pencil H — XJ in question is uniquely determined by the sequences 
{zfcj^o and Wk '■= fizk), k e Z+. So, as we already mentioned, the following question 
naturally arises. 
Nevanlinna-Pick problem. Is the function (f e Rq satisfying the interpolation relation 

(7.1) ip{zk) = Wk, keZ+ 

uniquely determined by the data {^fcl^g, {wk}'kLo'^ 

More details about Nevanlinna-Pick problems can be found in |[T], 1271 . l36l . 

Remark 7.1. Recall that an R-function is a function which is holomorphic in the open 
upper half plane C^ and maps C+ onto C+. For convenience, it is supposed that every 
(^ G R is extended to the lower half plane C_ by the symmetry relation ip{X) = <^(A), 
A e C_. Clearly, Rq is a subclass of R. In fact, the condition ip G Rq means that ip is an 
R-function and satisfies the following tangential interpolation condition 

(7.2) ^^^^^'\^''{x)' ^^'^• 

Roughly speaking, ( 17. 2b can be interpreted as the interpolation conditions (p{oo) = 0, 
ip'{oo) = —1. So, the Nevanlinna-Pick problem in question is a sublass of Nevanlinna- 
Pick problems in R. 

Before answering the question of the Nevanlinna-Pick problem we will prove the fol- 
lowing auxiliary statement. 

Lemma 7.2. We have that for j e Z+ 

ea = {H- ZjJ){C[o,j] (zj) + m[oj] (zjOtt^o.j] {zj)) = 

= iH[o,j] - ^jJlo.]]){£.[Oj] i'zj) + m[o.j] (zi)7r[oj] (zj)) 

11 fii I'l^^ 

Moreover, ifJ~2HJ^2 is self-adjoint in £^ then the systems < {J^^HJ^^ — 2^)^^ J^seo > 

I J j=o 

and {J^ej}'j°^Q are equivalent, that is, 

span{{J^^HJ^^ — zo)~\ . . . , {J~^HJ^^ — Zfe)~^eo} = spanjj^eo, . . . , J^Ck} 
for every k G Z_|_. 
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Zj<^j 



0. Then it follows from ( |43] |. and i4A\ that 



Proof. Notice that bj 

(7.4) 

(i7-ZjJ)7r[oj](zj) = {H[o,j] -ZjJ[oj])7r[oj](zj) = -{bj -Zjdj)Pj+i{zj)ej, 

{H-ZjJ)^[a.]]izj) = {H[o,j] -ZjJ[o,j])^[o,j]izj) = -{bj-z-jd-j)Qj+i{zj)ej+eo. 
Now, (17.31 ) is immediate from (17.41 ) by taking into account 

Qj+T-{Zj) _ Qj+l{Zj) 



Hoj]izj) 



Pj + l{Zj) 



■^3+iyzj) P,+i(z,) 

If J^^HJ^"! is a self-adjoint operator in £^ then ( 17.31 ) implies that 

(7.5) {J^^HJ^^ - Zj)"V"2eo =; J^(^[oj](^j) + TO[o,i](^j)7i'[Oj](2j)). 



Now, the equivalence follows from ( 17.51 ) for j = 0, . . . , fc and the fact that Qj{zj) + 
r?T,[o j](zj)Pj(zj) 7^ for j ~ 0, . . . ,k. The latter fact immediately follows from ( |3.2| i. 
(|4.13t , and the Liouville-Ostrogradsky formula ( 14. 2t . D 

Proposition 7.3. /ff/ie operator J^^ H J^^ is self-adjoint in l^ then the corresponding 
Nevanlinna-Pick problem ( 17.1b has the unique solution 

(/9(A) = m(A) := {{J-^HJ-^ - A/)~V~^eo, J~5eo). 

Proof. Clearly, for every A e C+UC- there exists a sequence r„ (A) g span{ Jaeo, . . . , J2e„} C 
dom(j2if Ji) such that 



(7.6) ||(J-5HJ-2 _ AK(A) - J-5eo|j ^ 0, 

It follows from Lemma|T2]that 

(7.7) 



oo. 



,(A) = Y. CkiX)iJ-^HJ-i - Zk)-'jho. 



fe=0 



Further, let HJ ^ = /g tdEt be a spectral decomposition of J 2 iJ J 2 . Then the func- 
tion 



i{X) 



t — A 



is a solution of the Nevanlinna-Pick problem (17.1b . Really, according to ( 17.3b we have 

m(zj) = ((J~5iyj-2 - Zj)~V^5eo, J-2eo)f2 =TO[oj](zj). 

Further, due to ( 12. 7b and ( 14. 13b one easily gets that m{zj) = Wj for j E Z_|_. Suppose that 
there is another solution fp{X) = /^ -frj-' Then we have 



it 



' '^ t-z, 

k=0 J 



dp{t) 



fc=0 "' 



d{Eteo, eo) 



\{J-^HJ-i - X)J2ck{X){J-^HJ-i -^j)-^J-ho - J^^ 



eol: 



fc=0 
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as n — >■ oo. Now, l/(i — A) is bounded for i e M since A e C+ U C_. Thus 

2 

dp[t) — > 0, n — > oo. 



Cfc(A) 1 



^ t-zk t- X 

k=0 



Finally, it follows that 



ipp{X) = lim / > - — —dp{t) 
"^°° -^^t^o '' 

is independent of p. Since (pp determines p (see for instance H] Chapter III]), all p's must 
be the same. D 

Proposition 7.4. If the operator J^^i HJ^^ is not self-adjoint in (? then the correspond- 
ing Nevanlinna-Pick problem (17. It has an infinite number of solutions. 

Proof. Since the deficiency indices of J^^HJ^^ are equal it has self-adjoint extensions 
in £^. Let Hi and H2 be two different self-adjoint extensions of J^^ HJ~2 in £'^. Then 
the following two functions 

ifiiX) = {{Hi - X)-'j-ho, J-5eo), ip2{X) = {{H2 - A)-V-5eo, J-^eo) 
are solutions of (17.1b . Really, according to Lemma l772l we have 

fk{'Zj) = {{Hk -Zj)"V"5eo, J"3eo) = {[H^oj] - Zj J[o,j])""^eo, eo) =% 

for every j e Z+ and k = 1,2. Since <y9fe G Ro, one also has (pk{zj) = Wj. 

Further, let A S C+ \ {zj}°^^. Note, that go = J^^cq ^ ran(J"^i7J"^ - A). 
To see this, suppose the contrary that there exists x e dom(J^2i7J^2 — A) such that 
go == {J-iHJ-^ - X)x and that {{J-^HJ-^-)*-X)y = 0. Then 

(50,2/) = {{J~--HJ--^ - X)x,y) - (x, ((J-^i7J-^)* - A)?/) = 0. 

We thus see that {go,y) = and {{J^^HJ^^)* — X)y = 0. As a consequence, the 
coefficients Uk = {gk,y) of the vector y ~ Tl'kLo^'-kfk solve ( 13.1b with the initial con- 
ditions u_i = ug = 0. Therefore, y = 0, that is, J^^HJ'^ is self-adjoint in £^. By 
hypothesis, this is false, so J^^eo ^ ra.n{J^^ HJ^2 — A). Thus (iJi — A)^^ J^^eo 
and (i72 - A)^V^2eo are in dom((J^2i7J-2)*) \ dom(J^2iJJ-2). So, we have 
{Hi — A)^^ J^2eo 7^ (i?2 — A)^^ J^2eo because otherwise, according to the fact that 
J~^HJ^2 has deficiency indices (1,1) and the von Neumann formulas we would have 
Hi = H2. 

Letr/ = (iJi-A)-ij-5eo-(iJ2-A)"V-5eo. Then one has ((J-^iJJ-5)*-A)77 = 
and, so, the coefficients rfk = {gk, v) of the vector 77 ^ J^kLo Vkfk give a solution 
of (13.1b with the initial conditions 

77-1 =0, 770 = {go,v)- 

Since 77 7^ we get (go, v) 7^ 0. As a consequence, we have ipi ^ ^2- To complete the 
proof it remains to observe that the function 

<^q(A) = aipi{X) + (1 - a)ip2{X) 

is also a solution of (17.1b for every a £ (0, 1). D 
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Remark 7.5. It follows from the proof that every self-adjoint extension of the symmetric 
operator J^^ HJ^2 generates a solution of the corresponding Nevanlinna-Pick problem. 
Moreover, by using the standard technique of theory of extensions of symmetric operators 
(see im, 1241 . 1391 ). one can get the description of all solutions of the Nevanlinna-Pick 
problem and it will be done elsewhere. The description of all solutions can be found, for 
instance, in 



The following theorem immediately follows from Propositions 17.31 and 17 .41 

Theorem 7.6. The Nevanlinna-Pick problem (17. Il l has a unique solution iff the corre- 
sponding operator J^^HJ^^ is self-adjoint in (? . 

Remark 7.7. Other criteria for the Nevanlinna-Pick problems to be determinate can be 
found in ll27l . ll36l . It is worth noting that, in the matrix case, the Stieltjes type criteria for 
Nevanlinna-Pick problems to be completely indeterminate were obtained by Yu. M. Dyu- 
karev in his second doctorate thesis (see 125)1 . l26l ). 

8. Convergence of multipoint Pade approximants 

In this section we prove a Markov type result on convergence of multipoint diagonal 
Pade approximants for Rq -functions. 

At first, let us recall that for the symmetric matrix J,^ \ ^[o. j] 'Ag ^1 '■^^ following esti- 
mate holds true 

(8.1) ll(^M^[o..]'^J]-^)"'ll<^^' ^eZ+. 

Before showing the convergence result, it is natural to obtain the precompactness. 

Proposition 8.1. The family {n^[o.j] }?^o " precompact in the topology of locally uniform 
convergence in C \ M. 

Proof. Let us rewrite the function m[o j] as follows 

"^[ojl(A) = ((./[oJ]^[ojy[oJ] - ^)"''^[oJ]eo,^[oJ]eo)- 
It follows from the Cauchy-Swarz inequality and ( |3.8l l that 

(^[oj]eo,eo) 1 



< 



t which, in view of the Montel theorem, implies the precompactness of {?ti[o ^j }°^o- '-' 

Now we are ready to prove the main result of this section. 

Theorem 8.2. Let a sequence of distinct numbers {zj}^^ C C-|- be given and let tp be 
a unique solution of the Nevanlinna-Pick problem ( 17. II ). Then all the multipoint diago- 
nal Pade approximants for (p at {zq, zq, . . . , Zj, Zj, . . . } exist and converge to (p locally 
uniformly in C \ R. 

Proof. Proposition 14.131 savs that the rational function mjoj] is the (JH-l)th multipoint 
diagonal Pade approximant. Further, according to Theorem 17.61 one obviously has that 
J^2i/J^2 is self-adjoint in £^ and, therefore, (J^siJJ^s — A)^^ is bounded for A e 

C \ R. Let -0 be a finite sequence, that is, ■0 = (-01, . . . , f/)/;, 0, 0, . . . )^. Then 

(i/-AJK' = (if[oj]-AJ[o,j])^ = <^ 
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for sufficiently large j E Z+ and (p is also a finite sequence. Further, one obviously has 
(8.3) 
[{J-^HJ-^ - Xr\rh, J-'eo) = ^Inn (( J-j]H[o,,] J"! - A)-^;^]]^/', J[olfo) 

In particular, formula ( 18.3b is valid for 

where r„ is defined by ( 17.7b . So, due to ( 17.61 ) we have that 

(8.4) J^20^ ^ j-2eo as n ^ oo. 
Moreover, the vectors (/>„ satisfy the following relation 

(8.5) "^[oj]'?^" ~^ "^[oJl^o as n -> oo 
for j e Z+. To see the latter relation, note that ( 18.4b implies 

(J~2 0„^7y) ^ (j-2eo,77) as 71^00 

for every t] e ^^. Putting 77 = J"^ Jiq '^i^k, k = 0, ... ,j, we get ( 18.51 ) from the fact that, 
in finite-dimensional spaces, the weak convergence is equivalent to the strong one. Now, 
taking into account ( 18. lb . ( 18.3b . (18.4b . and ( |8.5b . we obtain that ( 18. 3b holds true for (f) = eg, 
that is, 

for any A e C \ M. Finally, the statement of the theorem follows from the precompactness 
and the VitaU theorem. D 

Remark 8.3. In the case when (p G R[q;, /3] and the interpolation points stay away from 
[a, /3], an analog of the Markov theorem for multipoint diagonal Fade approximants is 
well known 1301 . 1471 (see also l23l where the operator approach was presented). In the 
case when the interpolation points belong to [—00, 0), the locally uniform convergence of 
multipoint Fade approximants for (p G R[0, +00) was proved under the Carleman type 
condition l37l (see also l38l where results in this direction are reviewed). It should be also 
remarked that there are some results on convergence of multipoint Fade approximants for 
rational perturbations of the Cauchy transforms of some complex measures Q, JD. 

It is a standard fact that the following condition 

Imzfe 



implies the determinacy of the corresponding Nevanlinna-Pick problem in Rq ll27l . l36l . 
Thus, the underlying operator J^ 2 iJJ^ 2 is self-adjoint in £'^. 

Corollary 8.4. If the given sequence {zjj^^o ^'^ti^fi^'< ( 18.6b then for every ip e Rq all 
the multipoint diagonal Fade approximants for Lp at {zq, zq t ■ ■ i Zj,'Zj, . . .} exist and con- 
verge to Lp locally uniformly in C \ R. 



Remark 8.5. First, note that (18.6b is sufficient for the Nevanlinna-Fick problem in Rq to 
be determinate but not necessary (see ll27l Chapter IV, Example 4.2]). It should be also 
noted that, under the Szego condition and the negation of the Blashcke type condition, the 
locally uniform convergence of multipoint diagonal Fade approximants for (^ G R[a, /3] 
was proved in l46l (see also ||6l). 
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Now, we are also able to adapt Theorem l5.5l for the self-adjoint case. 

Theorem 8.6. If J^"! HJ^^ is self-adjoint in P' then for every A G C+ \ {^fcl^g there 
holds 



Y, |m(A)(Pfe(A) + dk-iPk-i{\)) + Qfe(A) + c)fe-iQfe-i(A)p = 



A;=0 



7Tl(A) — ?7l(A) 

A-A 



Proof According to Theorem 18.21 and ( 15.11 ). we have that Kj(A) -^ m(A) as j — > oo. 
Now, the statement directly follows from Corollary 15. 3l and the inequality ( |5.5l l. D 
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